Effect of trap anharmonicity on a free-oscillation atom interferometer 
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A free-oscillation interferometer uses atoms confined in a harmonic trap. Bragg scattering from 
an off-resonant laser is used to split an atomic wave function into two separated packets. After one 
or more oscillations in the trap, the wave packets are recombined by a second application of the 
Bragg laser to close the interferometer. Anharmonicity in the trap potential can lead to a phase 
shift in the interferometer output. In this paper, analytical expressions for the anharmonic phase are 
derived at leading order for perturbations of arbitrary power in the position coordinate. The phase 
generally depends on the initial position and velocity of the atom, which are themselves typically 
uncertain. This leads to degradation in the interferometer performance, and can be expected to limit 
the use of a cm-scale device to interaction times of about 0.1 s. Methods to improve performance 
are discussed. 

PACS numbers: 03.75.Dg, 37.25.+k 



I. INTRODUCTION 



Atom interferometry is a sensitive tool for metrology 
and probes of fundamental constants [HI2]- Traditionally 
it uses atoms in free space, but many groups have inves- 
tigated interferometry of trapped atoms, both thermal 
and Bose-condensed The use of confined atoms 

has several advantages, including the abilities to main- 
tain high density for interaction studies and to impose 
complex atomic trajectories via the trapping potential. 
Another significant benefit of confinement is that the 
atoms do not fall under gravity, so measurements can 
be extended to long times without requiring a long drop 
distance. 

One interferometer configuration that has been of 
recent interest is the "free-oscillation" interferometer 
(Fi g. [T ), in which atoms are confined in a harmonic trap 
[51 lllT ll3ffTo] . An off- resonant laser pulse is applied to 
the atoms and induces momentum kicks via Bragg scat- 
tering. Typically, the laser is a standing wave tailored to 
produce two wave packets with momentum kicks ±2hk 
for light with wavenumber k. The two packets are al- 
lowed to separate and complete a half or full oscillation 
in the trap, at which time the laser pulse is applied again. 
The fraction of atoms brought back to rest after the sec- 
ond pulse depends on the phase difference between the 
packets at the time of recombination, making the device 
an interferometer. 

The free-oscillation interferometer can be compared 
with the Michelson configuration [31 El [Zl IE] , in which 
the atoms are held in a potential that is nearly flat in one 
direction. The wave packets are separated in that direc- 
tion and one or more Bragg reflection pulses are used to 
reverse the atomic motion and bring the packets back 
together. The free-oscillation interferometer offers sev- 
eral advantages over the Michelson interferometer, chiefly 
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FIG. 1. Atomic trajectories x(t) in a free-oscillation interfer- 
ometer. Cold atoms are confined in a harmonic potential with 
period r. At time ( = 0a laser pulse splits the atomic wave 
function into two packets (solid and dashed curves) moving 
in opposite directions. After completing either a half or a full 
oscillation in the trap, the packets are recombined. (A full 
oscillation is shown.) The fraction of atoms that returns to 
rest depends on the phase difference developed between the 
two packets. 



stemming from the difficulty of achieving a perfectly flat 
potential [HE3 OH- 

However, it is not easy to achieve a perfectly harmonic 
potential either. Anharmonicity is likely to be an impor- 
tant limit to the performance of a free-oscillation inter- 
ferometer because it causes a phase difference between 
the two arms that depends on the position and veloc- 
ity of the initial unsplit packet. Thermal atoms will 
have a large spread in position and velocity, so anhar- 
monicity can be expected to severely limit the usability 
of a free-oscillation interferometer with non-condensate 
atoms. Even a Bose-Einstein condensate will have some 
variation in the initial conditions due to the uncertainty 
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principle, along with extraneous motional excitations 
from technical effects. 

In this paper, we calculate the phase shift of a free- 
oscillation interferometer that is induced by anharmonic- 
ity. We will ignore atomic interactions here, although 
they can in fact be important |15j . In Section Il| we cal- 
culate the anharmonic phase shift. In Section III] we es- 
timate typical magnitudes of the anharmonicity and sug- 
gest some techniques for mitigation of the effect. Finally, 
Section [TV] presents concluding remarks. 



II. PHASE CALCULATION 

A. Perturbative Approach 

With reference to the atomic trajectories shown in 
Fig. [IJ we take the atoms to travel along the x direction. 
We assume this is a principle axis of the trapping po- 
tential, and also that anharmonicity does not introduce 
significant coupling between x and the other directions. 
In this case, the relevant motion is governed by the one- 
dimensional Hamiltonian 
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where m is the atomic mass, p = mx is the momentum, 
ujq is the unperturbed oscillation frequency, and mf(x) 
is a perturbing potential, assumed to be small. 

The phase accumulated by a packet is determined us- 
ing the classical action [IJ [18] , 



m - 0(0) = 



i 



(T - V") dt 



(2) 



where the kinetic and potential energies T and V are 
calculated for the classical trajectory x c i ass (t). To be 
precise, this gives the phase of the packet wave func- 
tion ip(x,t) at the packet center (x) — x c i ass . We here 
ignore the possible impact of phase gradients across the 
packet, which require the inclusion of atomic interactions 
to correctly analyze [T31 E3 E] . 

The classical trajectory is determined using the equa- 
tion of motion 



ui n x = — 



dx 



(3) 



Perturbative techniques for the solution of ^ are well 
known [THl |20]. We write x — xq + x\, where the sub- 
script indicates the perturbative order. We take xq = 
Acos(ujt — 6), where w = ujq -\-uj\ includes an amplitude- 
dependent frequency shift uj\ . Substituting these expres- 
sions into ([3]) and collecting the first-order terms yields 



±\ + UIqXi = 2oJqU>iXo 



dx 



(4) 



The frequency shift uj\ is determined by requiring that 
the right-hand side of Q contains no terms oscillating 



at wo, since these would drive unbounded excitation of 
x\. To first-order, this yields 
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For the interferometer phase shift, we consider the case 
where atoms make one full oscillation in the potential 
before recombination, as shown in Fig. [I] The phase 
developed by a wave packet during this motion is, from 
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,2^2 



- /(*) dt 
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where r is the period. Using x — xq + x\ and neglecting 
terms higher than first-order yields 



m 



1 



(±1 ~ ufal) ~ f(x ) 



dt. 
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The terms involving XqXi vanish upon integration be- 
cause xq and x\ are, by construction, orthogonal on this 
interval. To first-order, the derivative term in |7]) can be 
expressed as 



Xq = (cJq + 2uj lo 1 )A 2 sm 2 (wi - 6). 



(8) 



Using this to evaluating the integral leaves a perturbation 
phase 
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f(x ) dt. 



(9) 



Two cases can now be considered. The first supposes 
that f{x) varies over a length scale that is small compared 
to the amplitude of the atomic motion A. This might 
be caused by speckle in a laser trap or roughness in the 
conductors of a magnetic chip trap. In this case, u>\ will 
be small because df /dx in Eq. ^ will rapidly oscillate. 
In this limit, the phase shift will be dominated by the 
second term in ([9|, which is simply the integral of the 
perturbing potential: 



h->—jrJ f(x a )dt. 



(10) 



Further analysis would require knowing f(x), which will 
in general be specific to a particular apparatus. 

The second case supposes that f(x) is slowly varying 
compared to the atomic amplitude. This might result 
when the confinement potential is created by a distant 
magnetic or optical element. Here it is reasonable to 
Taylor expand f(x) as a power series 



(ii) 



n— 3 



where R is a characteristic length scale, which can typi- 
cally be taken as the distance to the trapping element. If 
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i<i?, then the dominant term will be the lowest power 
of n for which f n is non- vanishing. In a symmetric trap, 
/„ will be suppressed for all odd n, and careful design 
may result in the suppression of /„ for one or more even 
n as well. In general, however, u>\ will be non-zero and 
both terms of Eq. ^ will contribute. This case corre- 
sponds to the conventional anharmonic oscillator, and it 
is the main focus of the present paper. We continue its 



analysis in Section II B below. 

In the intermediate case, where the length scale of / is 
comparable to A, the effect on the phase is complicated. 
If the detailed form of f(x) is known, then Eqs. ^ and 
([9]) can be used to numerically estimate the phase shift. 
However, direct numerical computation of Eq. ^ would 
require similar effort and give greater accuracy. 



equal to the periods for the two packets, but differs by 
small amounts St± — t c — 2tt/lj±. In terms of the fre- 
quency shifts u>x±, this implies 

2tt 



St, 



St- 



(Wl+ - UJi-) . 



(16) 



Since the St's are already first-order in A n , their effect 
on the trajectories and phases can be calculated in zc- 
roth order. The trajectory is x±(t) — x a cosuj±t + (u a ± 
Mo) sinw+i, so setting x+ — X- at t c yields (u a +uo)St + = 
(u a — Uo)St-, leading to 

ir\ n h n 



St+ = 



-^±1 \ (A n r 2 - A n _~ 2 ) . (17) 

u ) 

The zero-order phase shift developed over a short time 
St is readily calculated from the action to be 



B. Anharmonic Trap: First Order 

As explained above, we here consider the case of a 
power-law perturbation which we will express as 
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An X n 
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If n is odd, then the result is simple: u)\ in Eq. ^ and 
m Eq. (9) are both zero by symmetry. There is thus 
no first-order effect on the interferometer, but we will 
consider the second order effect in Section Hi CI 

If n is even, then the frequency shift can be expressed 

as 



oil = 
where 



— A n - 2 f cos n (wt -6)dt= —A n ~ 2 h r 
2tt J cj 
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Using this in Eq. ([£]) yields a trajectory phase 
= I 1 - — ) '^-A n \h n (n even). 



(13) 
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Equation (151 gives the first-order phase acquired by a 



wave packet during one complete orbit through the con- 
fining potential. However, in the actual interferometer, 
neither packet makes a perfectly complete orbit. Sup- 
pose the wave packet starts at mean position x a and 
mean velocity v a . It will be convenient to work with 
the scaled velocity u a = v a /uJo- The initial packet is 
split into two, with one packet acquiring a velocity im- 
pulse +uq = 2hk/(mu>o) and the other —uq. The two 
packets therefore have motional amplitudes A± given by 
A\ = x 2 l +u 2 l +u 2 ] ±2u a uo. Because these amplitudes are 
different, the packets will experience different motional 
frequencies u± = luq + w\±. The interferometer will be 
complete when the two packets cross at time t c , with 
x+(£ c ) = X-(t c ). The crossing time is approximately 
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-A 2 St cos 20. 



(18) 



Here A\ cos 26 ± = x 2 — (u a ±u ) 2 . The additional phase 
difference from the trajectories is then 

- ^- - -^^(^r 2 - ^- 2 ){x 2 a +ui- uD 

(19) 



This adds to the the result ( 15 ) to give the total anhar- 



monic phase for a full-cycle interferometer, 

01+ + ^0+ — $1- ~ ^ < P- ( n even ) 
irm\ n h n 



1 - 



;-) (A n + -Ai) 



A r : 



(20) 

We note that in practice, the determination of t c is sub- 
ject to experimental uncertainty, which can lead to ad- 
ditional phase shifts. This issue is discussed further in 
Section HHC] below. 

From Eq. (20), the first-order anharmonic effect can 



be calculated for any power n. In particular, it is evident 
that a nonzero phase shift can be obtained only when 
v a = ujQU a is nonzero, since otherwise the amplitudes A + 
and A- will be equal and both packets will trace out the 
same orbit through the trap. The same argument can be 
made for the more general result of Eq. Q . 

The values of n most likely to be of interest are 4 and 
6. Evaluating (20) for these cases yields 
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u u a (x a +u a + u ). 



(21) 



(22) 



We checked these calculations by comparing to numer- 
ical integration of the action (|6| using the exact equation 
of motion (|3| . Figure [2] shows a comparison of the nu- 
merical and analytical results. In all cases examined, we 
found good agreement. 
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For later use, we note that the same result holds for even 
n, except that u>q — > u> and the sum is over odd s > 1. 

Extending the perturbation series to second order, we 
have x = xq + xi + X2 with u) = loq + uii- Inserting into 
the equation of motion Q and collecting second order 
terms gives 



x 2 + ujqX 2 = -(n - l)XxiA n 2 cos" 2 (uj t - 
+ 2uiQU>2^cos(ujQt — 9). 



(27) 



Again, the right-hand side must have no loq components. 
After some manipulation, this leads to 



U>2 = 



(n-l)A^L 2 



2^ 3 



E 



i 



(l + 6 s0 ). (28) 



FIG. 2. Comparison of first-order analysis (dashed curves) 
to numerical calculation (solid curves) of the anharmonic 
phase shifts in a free-oscill atio n interferometer. Here A<f> = 
(a>o/wr)A0, with A<f> from (20l. The recoil frequency for the 
Bragg laser beam is uj r = mvo/8h. The dimensionless an- 
harmonic coefficient A„ is given by X„Vq~ 2 /loq . The curves 
are labeled with the power of the anharmonic perturbation n. 
Both cases shown use v a = O.Ivq and x a = 0.1v /uj . 



To calculate the phase shift, we expand the Lagrangian 
to second order and discard terms that trivially integrate 
to zero. This yields 
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i / • 2 2 2\ \ n— 1 
+ -{x 1 -uJ x 1 )-X n x Xl 



dt. 
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C. Anharmonic Trap: Second-Order 

As noted previously, there is no first-order phase shift 
for a power-law perturbation with odd power n. How- 
ever, if an odd-n coefficient A„ is relatively large, then 
the second-order effect from it could be comparable to 
the first-order effect from an even power at higher n. In 
particular, terms scaling with (A3) 2 have the same am- 
plitude dependence as those scaling with A4, so a second- 
order calculation of the odd-n terms is needed to make a 
fair comparison of their impact. 

To proceed, we shall require the order-one correction 
to the motion, X\. Its equation of motion is 

fi+wgaji = -KA 71 - 1 cos"" 1 (wot- 0). (23) 

The cosine term here can be expanded using 



{cos u) p = Y, h< f ] 

cos su, 



(24) 



where the sum is over odd s from 1 to p when p is odd, 
and over even s from to p when p is even. The expansion 
coefficients are 



pi 



(p + s)!!(p-«)!! 



(2 - S s0 ). 



(25) 
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(123} is then 



Here h is identical to h n in Eq. ( 14 ) . The solution to 

1 cos s(uj n t — 9). (26) 



XnA n-l g h in-D 



1 1 / ■> -1 

Wn ^ S Z — 1 

u even s 



The first pair of terms are evaluated just as for even n 
to give a contribution 4>a — TrmA 2 u)2/Ti- The second 
pair of terms can be evaluated using expression ( 26 ) and 
orthogonality to obtain 



7T m\ 2 A 2n - 2 

7T mA 2 Ul2 

n-l h ' 



^ [h [ r 1] ] 2 

1^ l_ s 2 



(1 + S s0 ) 



(30) 



The final term can be evaluated similarly by expanding 
the cos" _1 (woi — 9) factor to yield 



27T mA 2 w 2 



•) C = 



(31) 



n — 1 Ti 

Combining these three terms gives the net phase for odd 



(„) 7r(n — 2) mA 2 uj2 



1 



(n odd). 



(32) 



with UJ2 from Eq. ( 28 ) 



To (32) must be added the correction Scf> accounting 



for the difference in oscillation frequencies for the two 
trajectories in the interferometer. This proceeds exactly 
as for even n. In fact, the leading-order phase shift for 
any n can be expressed as 



Acj) 



(n) _ 



— [c„ (5uj + A\ - 5lo-A 2 _) 

- ipl + u 2 a - ul) (Su + - Suj-) 



(33) 
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FIG. 3. Comparison of leading-order analysis (dashed) to 
numerical calculation (solid) of the n = 3 anharmonic phase. 
The dimensionless quantities A(f> and A3 are defined as in 
Fig. [5] with A</> from Eq. (341. The calculation uses v a = 
0.1«o and x a — O.luo/wo- 



FIG. 4. Comparison of second-order fit (dashed curves) to 
numerical calculation (solid curves) for a full-cycle interfer- 
ometer with quartic (n = 4) anharmonicity. Plotted values 
are scaled as in Fig. [2| with A4> from Eqs. (2l) and Q using 
v a = O.luo and x a = 0.1«o/wo- 



where for even n, c n = (n — 2)/n and 8lo ~ while 
for odd n, c n = (n — 2)/(n — 1) and Suj — u; 2 . The most 
interesting odd case n = 3 gives 



A<4 3) 



107r mA§ 3 



(34) 



which is similar to Eq. ( 21 ) as claimed. We again checked 



this result against direct numerical calculation of the 
phase and found good agreement, as seen in Fig. [3] 

In principle, the analytical approach used here could 
be continued to provide the next-leading-order correction 
to the phase. However, the calculation rapidly becomes 
complicated (2TJ [22] . We did, however, numerically in- 
vestigate the second-order correction for the quartic po- 
tential, as this seems to be the case most likely to have 
practical importance. We find a net interferometer phase 
of 



A# = -K 



KiU 2 a 



K 2 X 2 a 



K 3 X a U a 



(35) 

with dimensionless K = 20.61, k± = 1.667, k 2 — 1.001, 
and k 3 — 2.157. The estimated accuracy for A</> 2 is 0.2%, 
based on the consistency of the fit across a range of pa- 
rameter values. Figure [4] shows an example comparison 
including both first and second order contributions. We 
note that ( 22 ) and ( [35] ) have similar scaling with ampli- 
tude, so if neither A4 nor Xq are suppressed, then both 
terms must be considered if this level of accuracy is re- 
quired. 



also be implemented using just one half oscillation, with 
the packets recombined at t s» t/2. Here the phase will 
clearly be more sensitive to asymmetries in the confining 
potential, but it also becomes sensitive to asymmetric 
effects that might be of interest. It is worth noting, how- 
ever, that a half-cycle interferometer remains insensitive 
to uniform static forces such as gravity, since a constant 
force will simply shift the center of the harmonic poten- 
tial with no effect on the trajectories or phase. 

For even n, the phase calculation proceeds just as in 
the full-oscillation case, but all terms are reduced in mag- 
nitude by a factor of two. The final phase shift is there- 
fore just one half of Eq. (20 1. 



The effect of odd anharmonicities changes more signif- 
icantly, as expected. There is now a first-order contribu- 
tion to the phase, 



An) 



T/2 X 



dt 



that is readily evaluated to 



(„) _ 2mA w („) 



(36) 



(37) 



To this order, there is no correction to the oscillation 

(n) 

frequency, so <p\ la {f is the only contribution to the inter- 
ferometer phase difference. Here A±shi9± — u a ± u . 
Evaluating 0half for the case n = 3 yields 



D. Half-Cycle Interferometer 



2mA 3 , 2 j ., 

uo {u + x a + 2,u a ) 



(38) 



Until now we have supposed the atoms to complete a 
(nearly) full oscillation in the trap. An interferometer can 



Here the lack of symmetry means that a non-zero effect 
is obtained even for u a = 0. 
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III. DISCUSSION AND IMPACT 

A. Estimation of Effect 

Evaluating the anharmonic phase shift for a particu- 
lar experiment obviously requires knowing the value of 
the lowest-power unsuppressed A„ coefficient, either by 
measurement or by calculation from the trap geometry. 
Typically, however, an order-of-magnitude estimate for 
A ra can be obtained from the power-series expansion of 
Eq. (11). For the appropriate length scale R, the co- 



efficients /„ in that expansion should have comparable 
magnitudes. Furthermore, in the typical case that the 
harmonic confinement is provided by the same elements 
that introduce f(x), the expansion could be extended to 
n = 2 and the ji coefficient, given by fi = ujqR 2 /2, will 
have comparable magnitude to the other / n 's. In terms 
of the A„ coefficients of Eq. (12), this implies 



A,, 



2R n - 



(39) 



This allows an estimation as to when anharmonic ef- 
fects are likely to be important. For instance, a symmet- 
ric trap (with A dd ~ 0) would have 



(4) 



6tt 



mv a v% 
huj$R 2 



(40) 



Using Rb atoms with u>q = 2tt x 10 Hz, vo — 1.2 cm/s, 
and R = 1 cm, maintaining Acj) -C 1 requires the ini- 
tial atomic velocity v a <C 0.5 mm/s. The same con- 
figuration with w reduced to 2n x 1 Hz would require 
v a <C 0.5 //m/s. The strong dependence on uj reflects the 
fact that a weaker trap will allow the atomic trajectories 
to extend to larger distances where the anharmonicity is 
more significant. 

In principle, a fixed phase shift could be measured and 
subtracted out, but in practice v a is likely to fluctuate, 
making the anharmonicity into a source of noise. For 
thermal atoms, the velocity spread will be determined 
by the gas temperature, while for condensate atoms it is 
limited by the uncertainty principle and technical effects. 
In our experiments (23) . we produce Rb condensates in 
a relatively tight trap and then adiabatically reduce the 
confinement to give an oscillation frequency in the range 
of 1 to 10 Hz. After this process, we typically observe 
a center-of-mass motional excitation corresponding to a 
velocity variance cr 2 of about u>o x 10 -8 m 2 /s. We at- 
tribute this to a combination of imperfect adiabaticity, 
forces from uncontrolled ambient magnetic fields, and 
mechanical vibrations of the apparatus. In comparison, 
the non-interacting harmonic oscillator ground state has 
ct 2 /w = h/m = 7.3 x 10 -10 m 2 /s. The fundamental 
velocity uncertainty can be even lower in an interacting 
condensate |24j , but the interplay between anharmonicity 
and interactions requires additional consideration beyond 
the scope of this work. 



In practice, then, anharmonic effects can be expected 
to limit the usable oscillation frequency for a given trap 
geometry, and thus the measurement time r of the inter- 
ferometer. In the case of quartic anharmonicity with our 
empirical a v , the phase fluctuations will reach one radian 
at r « (0.09 s)i? 4 / 5 , for R in cm. This increases by about 
a factor of two for the ideal ground state a v . Anharmonic 
effects can thus be expected to impact a cm-scale device 
operating with interaction times greater than about 0.1 s. 



B. Amelioration 

An obvious way to reduce the impact of trap anhar- 
monicity is to reduce the anharmonicity itself by using 
a larger trap geometry to increase R. However, practi- 
cal applications often favor a more compact apparatus. 
Also, in both magnetic and optical traps the electrical 
or optical power required increases rapidly with the trap 
size. 

Alternatively, a small trap could be designed with re- 
duced anharmonicity. Typically this would be achieved 
by tuning one or more A„ coefficients to zero. For ex- 
ample, with our empirical a v a trap with A4 ps and 
A 6 ps 3w 2 /i? 4 would give one radian of phase noise at 
r ps (0.5 s)i? 8 / 9 , again with R in cm. This is about five 
times better than that obtained with unsuppressed A4. 
Controlling many coefficients in this manner, however, is 
likely to be challenging. 

We propose here another possibility for controlling an- 
harmonic effects: the use of a dual interferometer, as il- 
lustrated in Fig.[5j Here a wave packet initially at (x a , v a ) 
is split and allowed to propagate for one quarter oscilla- 
tion, after which the packets will be nearly at rest, with 
residual velocities v ps — uj x a for both. These packets 
are used as the sources for two independent interferom- 
eters: the packets are split, allowed to oscillate for one 
period, and then recombincd. The output signal ^ is 
taken as the difference between the phases of two inter- 
ferometers. Because the initial velocities of the two inter- 
ferometers are correlated, the leading-order anharmonic 
phases cancel in the difference. The phase difference is 
also less sensitive to mechanical vibrations and several 
other technical effects [25]. However, this configuration 
evidently requires the phase shift of mctrological interest 
to be differential between the two interferometers. For 
example, a local field might be applied to just a single 
packet. 

We analyze the dual interferometer for the n = 4 case. 
The trajectories are, from Eq. 



(26), 



x{t) = Acos(u}t — d) + — ^A 3 cos3(u;i - 9). (41) 
32l0q 

We use this solution to determine the actual initial condi- 
tions for the subsequent pair of interferometers. We here 
use + (— ) to label the interferometer derived from the 
packet originally given a positive (negative) momentum 
kick. 
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FIG. 5. Dual interferometer configuration. At time t = 0, 
a Bragg splitting pulse is applied to the nominally stationary 
initial wave packet in a nominally harmonic trap of period 
r. The resulting two packets move in the potential until they 
come to nominal rest at time r/4, when the splitting pulse 
is again applied. This produces two pairs of packets. All 
four packets are allowed to propagate for a measurement time 
r and then recombined with a final pulse. The upper pair 
of trajectories (solid curves) forms one interferometer with 
output phase A<f> + , and the lower pair (dashed curves) forms 
another with output phase A(j>-- The phase difference VP = 
A(j>+ — A<f>- features reduced sensitivity to the anharmonicity 
of the trap potential. 



The amplitudes and angles in (41 ) must be determined 



in terms of x±(0) = x a and u±(W) — u a ±uo. We write 
^4 = ^0 + ^! and 9 = 8q + 9\. In zeroth order, we 
have x a = Aq± cos 6q± and u a ± uq — Aq± sin 9q± ■ The 
first-order corrections are then calculated to be 



Ai-_ 



1 A 4 Ag ± 
32 w 2 
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COS Uq± 



9) 



(42) 



and 



-\^Al ± cos6 Q± sm.d Q± (3 + 2cos 2 (9 0± ) . (43) 



The frequencies ui± are given by ([13 



3 ^4 , 2 



1 o±- 



(44) 



After the first splitting pulse, the packets propagate for 
a time — 7r/(2w), where Co = (lu+ + w_)/2. Inserting 
this into the above trajectories yields 

u±(h) = -x a - —~ [37TM a u (u ± u a ) (45) 
8u 

+2x a (3uq + 3u 2 + x\ ± 6u u a )] . 

Using these as the initial velocities for the two subsequent 
interferometers in (21 1 gives a phase difference 



* = A(j> + - A0_ = - 



9ir uiXa 



4„,4 



4 TiojI 



u u a (nu a +4:X a ), (46) 



which is much smaller than the individual A</>'s. The 
result is second-order in A4, but the direct second-order 



effect of ( 35 ) also largely cancels to give a third-order 



correction to ^. The leading-order A6 effect (22) cancels 
as well. Using the estimated value A4 « 2ujfi/R 2 , our 
empirical values for a v , and taking a x a v /ujQ, we obtain 
a measurement time limit r ss 0.9R s/cm for < 1. 



C. Measurement of Anharmonicity 

Although the estimates above may be helpful, an ac- 
curate consideration of the effects discussed here will re- 
quire actual knowledge of the relevant anharmonic coef- 
ficients. While they can be calculated in principle, an 
experimental technique to measure them would likely 
be useful. The most straightforward approach is to 
observe how the oscillation frequency depends on the 
motional amplitude, via Suj. However, it may be dif- 
ficult to measure Su> with sufficient accuracy. For ex- 
ample, in a 10-Hz trap with quartic anharmonicity and 
R = 1 cm (A4 = 8 x 10 3 cm~ 2 s~ 2 ), Rb atoms with 
v a = 0.5 mm/s would experience a significant phase 
A(f> ss 1 rad. However, the frequency shift for atoms with 
amplitude A = vq/ujq w 200 /zm would be only 3 mHz. 
Such a small shift could be difficult to measure, given a 
finite lifetime of atoms in the trap. It may prove more ef- 
fective to use the interferometer to characterize the trap, 
by for instance measuring how the interference phase A<p 
varies with v a . The results presented here should be use- 
ful for this purpose as well. 

It should also be noted that imprecise knowledge of 
the oscillation frequency u can itself lead to phase uncer- 



tainty. In deriving ( 20 1 , we assumed a measurement time 



t c such that the two packet centers exactly crossed at the 
time of recombination. If the actual measurement time 
is too far from t c , then the interference contrast will be 
reduced because the packets will not be well-overlapped. 
For a small timing error, however, the overlap will remain 
large and the dominant effect will be a phase shift result- 
ing from the differing velocities of the two packets. This 
can be calculated from Eq. ( 18 ) along with a correction 



mv a Sx/h for packet separation Sx. The result is 

TTL 

5cj) = 2-v v a 5t (47) 



for timing error 5t — t — t c . Using our empirical value for 
the velocity uncertainty in Rb, this results in a significant 
phase uncertainty of (7 rad/ms)5t. As above, determin- 
ing uj with sufficient accuracy to avoid this problem may 
be challenging. This phase error can be reduced using 
the dual interferometer scheme, because to lowest order 
it is the same for both pairs of trajectories. 



IV. CONCLUSIONS 

We have calculated the leading-order effects of trap 
anharmonicity on a free-oscillation atom interferometer. 
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For a typical cm-scale device, anharmonic phase shifts are 
likely to be important for interaction times of about 0.1 s 
or greater, with the effects growing rapidly as the interac- 
tion time is increased. Possible methods for amelioration 
include nulling the low-order anharmonic coefficients via 
careful trap design, using a larger-scale trap with less 
anharmonicity, minimization of the initial velocity of the 



atoms, or phase cancellation in a dual interferometer. 

ACKNOWLEDGMENTS 

This work was supported by the National Science 
Foundation (Grant No. PHY- 0244871). 



P. R. Berman, ed., Atom Interferometry (Academic 
Press, San Diego, 1997). 

A. D. Cronin, J. Schmiedmayer, and D. E. Pritchard, 
Rev. Mod. Phys. 81, 1051 (2009). 

Y. Shin, M. Saba, T. A. Pasquini, W. Ketterle, D. E. 
Pritchard, and A. E. Leanhardt, Phys. Rev. Lett. 92, 
050405 (2004). 

Y. J. Wang, D. Z. Anderson, V. M. Bright, E. A. Cornell, 
Q. Diot, T. Kishimoto, M. Prentiss, R. A. Saravanan, 
S. R. Segal, and S. Wu, Phys. Rev. Lett. 94, 090405 
(2005). 

S. Wu, E. J. Su, and M. Prentiss, Euro. Phys. J. D 35, 
111 (2005). 

T. Schumm, S. Hofferberth, L. M. Andersson, S. Wil- 

dermuth, S. Groth, I. Bar- Joseph, J. Schmiedmayer, and 

P. Kriiger, Nature Phys. 1, 57 (2005). 

O. Garcia, B. Deissler, K. J. Hughes, J. M. Reeves, and 

C. A. Sackett, Phys. Rev. A 74, 031601(R) (2006). 

A. S. Arnold, C. S. Garvie, and E. Riis, Phys. Rev. A 

73, 041606(R) (2006). 

M. Horikoshi and K. Nakagawa, Phys. Rev. Lett. 99, 
180401 (2007). 

R. E. Sapiro, R. Zhang, and G. Raithel, Phys. Rev A 79, 
043630 (2009). 

P. M. Baker, J. A. Stickney, M. B. Squires, J. A. Scoville, 

E. J. Carlson, W. R. Buchwald, and S. M. Miller, Phys. 
Rev. A 80, 063615 (2009). 

F. Baumgartner, R. J. Sewell, S. Eriksson, I. Llorente- 
Garcia, J. Dingjan, J. P. Cotter, and E. A. Hinds, Phys. 



[13 

[14 
[15 
[16 
[17 
[18 

[19 
[20 
[21 
[22 
[23 

[24 

[25 



Rev. Lett. 105, 243003 (2010). 

J. H. T. Burke, B. Deissler, K. J. Hughes, and C. A. 

Sackett, Phys. Rev. A. 78, 023619 (2008). 

S. R. Segal, Q. Diot, E. A. Cornell, A. A. Zozulya, and 

D. Z. Anderson, Phys. Rev. A 81, 053601 (2010). 

R. P. Kafle, D. Z. Anderson, and A. A. Zozulya, Phys. 

Rev. A 84, 033639 (2011). 

J. A. Stickney, R. P. Kafle, D. Z. Anderson, and A. A. 

Zozulya, Phys. Rev. A 77, 043604 (2008). 

M. Olshanii and V. Dunjko (2005), arXiv: |cond-| 

mat/0505358l 

P. Storey and C. Cohen- Tannoudji, J. Phys. II France 4, 
1999 (1994). 

L. D. Landau and E. M. Lifschitz, Mechanics (Pergam- 
mon Press, Oxford, 1976), 3rd ed. 

V. Barger and M. Olsson, Classical mechanics: a modern 
perspective (McGraw-Hill, New York, 1995), 2nd ed. 
C. R. Eminhizer, R. H. G. Helleman, and E. W. Montroll, 
J. Math. Phys. 17, 121 (1976). 

A. I. Kuznetsov, O. V. Karagioz, and V. P. Izmailov, 
Meas. Tech. 48, 848 (2005). 

J. M. Reeves, O. Garcia, B. Deissler, K. L. Baranowski, 
K. J. Hughes, and C. A. Sackett, Phys. Rev. A 72, 
051605(R) (2005). 

J. Stenger, S. Inouye, A. P. Chikkatur, D. M. Stamper- 
Kurn, D. E. Pritchard, and W. Ketterle, Phys. Rev. Lett. 
82, 4569 (1999). 

G. T. Foster, J. B. Fixler, J. M. McGuirk, and M. A. 
Kasevich, Opt. Lett. 27, 951 (2002). 



